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We review the production and flow of identified hadrons at RHIC with 
a main emphasis on the intermediate transverse momentum region( px ~ 
2-5 GeV/c). The goal is to unravel the dynamics of baryon production and 
resolve the anomalously large baryon yields and elliptic flow observed in 
the experiments. 

1. Introduction 

This paper explores the relativistic heavy ion collisions at RHIC and the 
medium produced in these collisions using hadronic observables. Being most 
abundantly produced, hadrons define the bulk medium behavior, which is 
governed by soft, non-perturbative particle production. Analysis of identi- 
fied hadron spectra and yield ratios allows determination of the kinetic and 
chemical properties of the system. Hydrodynamics models |l] have been 
successful in reproducing identified hadron spectral shapes and their char- 
acteristic mass dependence at low pt&s well as the azimuthal anisotropy 
of particle emission. Notably, in order to match the data the models re- 
quire rapid equilibration of the produced matter and a QGP equation of 
state. The particle abundances also point to an equilibrated system and are 
well described by statistical thermal models [2J. The chemical freeze-out at 
T c h ~ 170 MeV is suggestive, as it is at the phase boundary of the transition 
between hadron gas and QGP, as predicted by lattice QCD calculations [3]. 

Above px ~ 2 GeV/c, hard-scattering processes become increasingly 
important. After the hard-scattering, a colored object (the hard-scattered 
quark or gluon) traverses the medium produced in the collision and in- 
teracts strongly. As a result, it loses energy via induced gluon radiation. 
This phenomenon, known as jet-quenching, manifests itself as suppression 
in the yields of high-py hadrons, when compared to the production in pp 
collisions and weakening of the back-to-back angular correlations between 
the jet fragments. The yield suppression is measured in terms of the nu- 
clear modification factor Raa = YieldAA/N co ii/Yield pp , where the number 
of binary nucleon-nucleon collisions, N co u, is introduced to account for the 
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nuclear geometry. In this paper, we use the ratio Rcp, which is obtained 
from the N co u scaled central to peripheral spectra and carries similar infor- 
mation. Jet quenching was discovered at RHIC both in suppressed hadron 
production at high-py [I] {Raa < 1) an d in vanishing back-to-back jet 
correlations |SJ. 

Another discovery, unpredicted by theory, is a large enhancement in the 
production of baryons and anti-baryons at intermediate px ~ 2-5 GeV/c |HJ 
[7j, compared to expectations from jet fragmentation. This is in contrast 
to the suppression of 7r° |2|. In central Au + Au collisions the ratio p/tt 
is of the order 1 - a factor of 3 above the ratio measured in peripheral 
reactions or in pp collisions. In this region of pt fragmentation dominates 
the particle production in pp collisions. It is expected that fragmentation 
is independent of the colliding system - hence the large baryon fraction 
observed at RHIC comes as a surprise. At RHIC, the medium influences 
the dynamics of hadronization resulting in enhanced baryon production but 
the exact mechanism is not yet completely understood. This paper reviews 
the latest experimental results relevant to this subject. 

2. Radial flow at intermediate px- 

The most common conjecture that is invoked to explain the large p/ir 
ratios observed by PHENIX 7_ is the strong radial flow that boosts the 
momentum spectra of heavier particles to high px- In this scenario, the 
soft processes dominate the production of (anti)protons at 2-4.5 GeV/c, 
while the pions are primarily produced by fragmentation of hard-scattered 
partons. In Fig. ^ we compare the 10% central spectra of tt^, K^,p, and 
p to a hydrodynamics model |U] that has been fitted to the data. The 
free parameters in the model are the kinetic freeze-out temperature Tj , 
the transverse flow velocity fir and the absolute normalization. The line 
drawn through the </>-meson spectrum is the model's prediction obtained 
after fitting all other particle species. We see that: 1) Hydrodynamics gives 
a good description of the p and p spectral shapes up to 3 GeV/c, and 2) 
the </>-meson spectrum can be described by the same parameter set as the 
protons. For lighter particles, the deviation from hydrodynamics happens at 
lower pt- These results may lead to the conclusion that the enhanced p/ir 
ratio is a mass effect and the intermediate p-r(a n ti) protons are primarily 
produced in soft processes. We now examine the scaling of the yields in 
different centrality classes. We expect that for soft production the yields 
will scale as the number of nucleons participating in the collision, while 
for hard processes the scaling is with N co u . In Fig. |2] the ^distributions 
for (anti)protons and are scaled down each by their respective N co u. To 
isolate mass effects from baryon/meson effects we compare a heavy meson 
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Fig. 1. Transverse momentum spectra of tt , K , and p,p and hydrodynamical 
fit results for 0-10% central Au+Au collisions at ^/sWn = 200 GeV [ID]. The 
Arranges for the fit are indicated by the solid lines, while the dashed lines show 
the extrapolated predictions for each particle species. The 0-meson spectrum, not 
included in the fit, is compared to the model's prediction. 



to the protons. The result is rather surprising. At intermediate pxthe p + p 
yields scale with N co u as expected for hard processes. The (p yields do not 
scale. Although the shape of the p + p and 4> spectra is the same and is well 
reproduced by hydrodynamics, the absolute yields for <p grow slower with 
centrality. When the central and peripheral yields are used to evaluate the 
nuclear modification factor (Fig. |SJ), the (anti)protons show no suppression 
(Rep ~ 1) , while the <fr are suppressed similar to ir°. This result rules 
out the radial flow (and the mass) as the sole factor that is responsible for 
the baryon enhancement. The similarity in the centrality dependence of <j> 
and 7T production suggests an effect related to the number of constituent 
quarks rather than the mass. The STAR experiment also observed a clear 
baryon/meson distinction in Rcp of K*, K®, A, and 3 |12j . 

3. Recombination and empirical scaling of elliptic flow. 

Recently, several quark recombination models |131ll4llTo] have been pro- 
posed to resolve the RHIC baryon puzzle. In the dense medium produced 
in central Au + Au collisions, recombination of quarks becomes a natural 
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Fig. 2. iV co /;-scaled transverse momentum spectra of (p + p)/2 and </>-mesons for 
three different centrality classes. High-p^ baryon yields scale with N co u, while no 
scaling is observed for 0- mesons |1U| . 



hadronization mechanism. When produced from the same underlying ther- 
mal quark distribution, baryons get pushed to higher prthan mesons due 
to the addition of quark momenta. At intermediate ^recombination wins 
over fragmentation for baryons, while mesons are still dominated by frag- 
mentation. After fitting the inclusive hadron spectra to extract the thermal 
component, the models are able to reproduce a large amount of data on 
identified particle spectra, particle ratios and nuclear modification factors. 
The most spectacular success of the recombination models comes from the 
comparison with the data on elliptic flow. At low-pj-hydrodynamics de- 
scribes both the magnitude and the mass dependence of t>2- However, at 
Pt> 2 GeV/c the mass ordering of V2 changes,namely : V2(p) > ^2(71") [TB] 
and i>2(A) > V2(K S ) [T12j . In addition, the size of the signal is too big to 
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be explained by asymmetric jet absorption ^Jj.The recombination models 
solve the problem by assigning the elliptic flow signal to the quarks, instead 
of the hadrons. Then the baryon/meson split in vi is naturally explained. 
It has been demonstrated empirically, that the flow per quark is universal. 
Recent results from the STAR ^H] experiment that include the measure- 
ment of multi-strange baryons are shown in Fig. ^ A clear baryon/meson 
difference is observed in the data at pr> 2GeV/c. The results from a typ- 
ical hydrodynamic model calculations ^5] are shown with a band. After 
re-scaling of both axes in Fig. 0]to represent the quark flow, the data falls 
on a universal curve as demonstrated in Fig. |SJ 



4. Jet correlations with leading baryons or mesons 

The recombination models resolve most of the baryon/meson effects ob- 
served in the data. However, from spectra, particle ratios and elliptic flow 
it is difficult to infer whether the recombining quarks come from the ther- 
mal bath (soft processes) or from hard-scattering. To unravel the nature 
of the baryon enhancement and to test the recombination approach, the 
PHENIX experiment examined the two-particle angular correlations with 
identified meson or baryon trigger particle |2Uj . The momentum of the trig- 
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Fig. 3. The nuclear modification factor, Rcp, of 0, p + p and ir° measured by 
PHENIX in Au+Au collisions at y/J^N = 200 GeV [TTl|T(?|. 
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gers was chosen to be in the range of pt where baryon/meson differences 
are observed (2.5 < pr < 4 GeV/c). For both types of trigger, clear jet- 
like angular correlations were observed both on the same side of the trigger 
and at 180°. This result shows that both mesons and baryons have a signifi- 
cant hard-scattering component at intermediate pr, although the p/ir ratios 
are dramatically different from fragmentation in the vacuum. In order to 
quantify the similarity between baryon and meson triggered correlations, 
the yield of associated particles is integrated at the near-side and the away 
side peaks. The results for different centrality classes and colliding systems 
are shown in Fig. 6. There is an increase in partner yields in mid-central 
Au+Au compared to the d+Au and p+p collisions. In Au+Au collisions, 
the near side yield per meson trigger remains constant as a function of cen- 
trality, whereas the near-side yield per baryon trigger decreases in the most 
central collisions as expected if a fraction of the baryons were produced by 
soft processes such as recombination of thermal quarks. The dashed line 
in Fig. 6 represents an upper limit to the centrality dependence of the jet 
partner yield from thermal recombination. 

The data clearly disagree with both the centrality dependence and the 
absolute yields of this estimation, indicating that the baryon excess has 
the same jet-like origin as the mesons, except perhaps in the highest cen- 
trality bin. The bottom panel of Fig. 6 shows the conditional yield of 
partners on the away side. It drops equally for both trigger baryons and 
mesons going from p+p and d+Au to central Au+Au, in agreement with 
the observed disappearance and/or broadening of the dijet azimuthal 
correlations. It further supports the conclusion that the baryons originate 



7 



0) 0.1 

O) 

O) 



2 

<D 0.05 



j_ -0.05 

<u 

.5> 0.1 



2 



0.05 



Near Side 



t f 



I 

! 



• meson triggers, AuAu 

■ baryon triggers, AuAu 

-| — i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i— | — i — i — i — i — i — i — h- 



Away Side 



o meson triggers, dAu 
n baryon triggers, dAu 
± charged hadron triggers, pp 



! 



1 



100 



200 



300 



N 



part 



Fig. 6. Yield per trigger for associated charged hadrons between 1.7 < pr < 2.5 
GeV/c for the near- (top) and away- (bottom) side jets [201 ■ The dashed line (top) 
represents an upper limit of the centrality dependence of the near-side partner yield 
from thermal recombination. 



from the same jet-like mechanism as mesons. The description of the data in 
the pQCD framework would require an in-medium modification of the jet 
fragmentation functions. For recombination models, the experimental re- 
sults imply that shower and thermal partons have to be treated on an equal 
basis [Tl]. 
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5. Summary 

We reviewed the results on hadron production and flow in relativistic 
heavy ion collisions at RHIC. The production mechanisms at low-p^and 
high-p^are relatively well understood in terms of soft and hard processes, 
respectively. The intermediate p^region (2 < px < 5GeV/c) is marked by 
a number of puzzling experimental observations and most notably, by the 
baryon excess over the expectation from vacuum fragmentation functions. 
By comparing spectra and centrality scaling of (anti)protons and 0-mesons, 
we established that the excess of anti-protons with respect to pions is not 
due to the larger mass of the anti-proton, but is related to the number of 
constituent quarks. The recombination models get a beautiful confirmation 
in the empirical scaling relation of the elliptic flow results. Jet-correlations 
with trigger baryons or mesons show a similar hard-scattering component in 
both. This observation is also in line with the N c oll scaling observed in the 
yields of protons and anti-protons. However, it implies that protons origi- 
nate from hadrons that experience little or no energy loss, while pions come 
from partons that have suffered large energy loss. This result is conceptually 
difficult, unless baryons and mesons have a very different formation time and 
thus - the original partons have different time to interact with the medium. 
Recombination models which combine hard-scattered partons with thermal 
ones give the most likely explanation of the experimental results as a whole. 
The baryon excess is clearly an effect of the medium produced in Au + Au 
collisions and maybe, thorough the comparison with recombination models, 
gives evidence for its partonic nature. 
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